The Spatial, Ionization, and Kinematic Conditions of the z = 1.39 
Damped Lya Absorber in Q0957 + 561 A,B^ 



Christopher W. Churchill^, Richard R. Mellon^, Jane C. Charlton^ 

Department of Astronomy and Astrophysics 
The Pennsylvania State University 
University Park, PA 16802 
cwc, rmellon, charlton@astro.psu.edu 

and 

Steven S. Vogt^ 

UCO/Lick Observatories 
Board of Studies in Astronomy and Astrophysics 
University of California, Santa Cruz, CA 96054 
vogt@ucolick. org 

ABSTRACT 

We examined the kinematics, ionization conditions, and physical size of the 
absorption clouds in a z = 1.3911 damped Lya absorber (DLA) in the double 
image lensed quasar Q 0957 + 561 A,B (separation 135hj^ pc at the absorber 
redshift). Using HIRES/Keck spectra (FWHM ~ 6.6 km s~^), we studied the 
Mgll AA2796,2803 doublet. Fell muhiplet, and Mgl A2853 transition in absorp- 
tion. Based upon the Fell profiles (the Mgll suffers from saturation), we defined 
six "clouds" in the system of sightline A and seven clouds in system of sightline B. 
An examination of the N{v) profiles, using the apparent optical depth method, 
reveals no clear physical connection between the clouds in A and those in B. 
The observed column density ratios of all clouds is log A(Mgl)/A(Fell) ~ —2 
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across the full ~ 300 km s~^ velocity range in both systems and also spatially 
(in both sightlines). This is a remarkable uniformity not seen in Lyman limit 
systems. The uniformity of the cloud properties suggests that the multiple 
clouds are not part of a "halo". Based upon photoionization modeling, using 
the iV(Mgl)/iV(FeIl) ratio in each cloud, we constrain the ionization parameters 
in the range —6.2 < logf/ < —5.1, where the range brackets known abundance 
ratio and dust depletion patterns. The inferred cloud properties are densities of 
2 < TT-H < 20 cm~^, and hne of sight sizes of 1 < D < 25 pc. The masses of the 
clouds in system A are 10 < M/Mq < 1000 and in system B are 1 < M/Mq < 60 
for spherical clouds. For planar clouds, the upper limits are AOOh^^ Mq and 
160h^^ Mq for A and B, respectively. We favor a model of the absorber in which 
the DLA region itself is a single cloud in this complex, which could be a parcel 
of gas in a galactic ISM. We cannot discern if the Hi in this DLA cloud is in a 
single, cold phase or in cold+warm phases. A spherical cloud of ~ 10 pc would 
be limited to one of the sightlines (A) and imply a covering factor less than 0.1 
for the DLA complex. We infer that the DLA cloud properties are consistent 
with those of lower density, cold clouds in the Galactic interstellar medium. 

Subject headings: quasars: absorption lines — quasars: individual (Q0957+561 
A,B) — galaxies: ISM 

1. Introduction 

Damped Lya absorbers (DLAs) are well appreciated as useful astronomical laboratories 
for measuring metal enrichment, nucleosynthetic processing, and dust evolution from z — A 
to the present epoch (e.g. Pettini et al. 1994; Lu ct al. 1996; Pcttini et al. 1997; Vlaldilo 1998; 
Pcttini et al. 1999; Prochaska & Wolfe 2000; Prochaska 2002; Ledoux, Bergeron, & Petitjean 
2002). Their redshift number density evolution (Wolfe et al. 1995; Storrie-Lombardi et al. 
1996; Rao & Turnshek 2000) places constraints on r2(Hl, z), and probes the Hi mass function 
for column densities greater than A^(Hl) = 2 x 10^° cm~^. The gas kinematics measured from 
low-ionization metal absorption lines in DLAs provides valuable data upon which hypotheses 
for their dynamical formation, evolution, and physical nature can be based. To date, these 
competing hypotheses include thick rotating disks (Prochaska & Wolfe 1997, 1998), merging 
cold dark matter halos (Haehnelt, Steinmetz, & Ranch 1998; McDonald & Miralda-Escude 
1999), and supernovae winds (Nulsen, Barcon, & Fabian 1998; Shaye 2001, however, see 
Bond et al. 2001). 

The fact that DLAs provide leverage on questions relating to cosmological chemical 
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evolution, nucleosynthetic and dust evolution, galaxy and/or hierarchical clustering evolu- 
tion, and global star formation evolution makes them one of the most useful, and therefore, 
important astronomical objects of study. Of the great deal that has been learned about 
DLAs, one property remains elusive; there is little direct observational evidence from which 
the characteristic size of the DLA "region" itself can be deduced. By DLA region, we mean 
the parcel of gas that is associated with the very large Hi column density greater than 
A^(Hl) = 2 X 10^° cm~^. Prom the size estimates, one could then deduce the covering factor 
of DLAs within their host objects, as well as their densities, and most importantly, their 
masses (though this quantity is geometry dependent). 

Associating DLAs with the normal, bright galaxy population and assuming a unity cov- 
ering factor, Steidel (1993) deduced a characteristic size^ of ~ 20h^^ kpc. However, it is now 
known that DLAs are associated with galaxies having a wide variety of morphological types, 
including 0.1 L* galaxies and low surface brightness (LSB) galaxies (Le Brun et al. 1997; 
Rao & Turnshek 1998; Bowen, Tripp, & Jenkins 2001). Indeed, the luminous host of many 
DLAs remain undetected even after dedicated searches to find their diffuse Lya emission 
(e.g., Deharveng, Bowyer, & Buat 1990; Lowenthal ct al. 1995), their optical continuum 
emission (e.g., Steidel et al. 1997; Kulkarni et al. 2000, 2001) and their Ha emission (e.g.. 
Bunker et al. 1999; Bouche et al. 2000; Kulkarni ct al. 2000, 2001). 

In the absence of a well-defined luminosity function for the luminous hosts of DLAs, 
and without knowledge of the covering factor of the DLA region in these objects, we remain 
ignorant of a characteristic DLA size. It is safe to say, however, that estimates incorporating 
low surface brightness galaxies in the luminosity function of known DLA hosts would drive 
down the characteristic size. Scaling the Steidel estimate by the ratio of for the luminosity 
functions of normal bright galaxies (Lilly et al. 1995; Ellis et al. 1996; Lin et al. 1999) and 
LSB galaxies (e.g., Dalcanton et al. 1997), we find an approximate characteristic size of 
(0.033/(0.033 + 0.08) • 20h^^ ~ Gh^^ kpc. 

The data directly constraining the physical sizes of DLAs are few. In the lensed quasar 
HE 1104 - 1805, Smette et al. (1995) placed hmits of 5-17h^^ kpc on a DLA in one of 
the images. A more constraining result was obtained by Michahtsianos et al. (1997), who 
used the lensed quasar Q 0957 + 561 A,B to show that the physical extent of a DLA at 
z — 1.3911 is less than 135hj^ pc. On the other hand, Petitjean et al. (2000) suggest a 
lower limit of 200/1^5^ pc for a probable DLA in the lensed quasar APM 08279 + 5255 (based 
upon partial covering arguments). Using VLBA observations of the slightly extended radio 
quasar B 0738 + 313, Lane, Briggs, & Smette (2000) find log A^(Hl) = 21.2 cm'^ in 21-cm 
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absorption at two locations separated by 20h^^ pc [we note, however, that 21-cm emission 
hne techniques may not be selecting objects that are one and the same as the DLAs selected 
using quasar absorption line techniques (see Rao & Turnshek 2000; Churchill 2001)]. 

In addition to the great insights garnished from low- to moderate-resolution spectro- 
scopic studies of multiply imaged quasars, high resolution spectra provide constraints on 

the kinematics of metal-line systems at the few kilometers per second level (e.g., Ranch, 
Sargent, & Barlow 1999; Lopez et al. 1999; Pctitjcan ct al. 2000; Ranch et al. 2002b). This 
means the gaseous structure can be studied component by component, directly providing 
their sightline to sightlinc velocity shears, At'/Ar, where r is the physical separation of the 
sightlines at the absorber. Such measurements allow insights into the ionization structures, 
masses, density gradients, and energy input budgets (e.g.. Ranch, Sargent, & Barlow 1999; 
Ranch et al. 2002a; DOdorico, Petitjean, & Christiani 2002). It is important to reahze that 
these insights are derived based upon models of the gas dynamics and geometry. These 
models, however, arc strongly founded upon observations of the gaseous components in the 
Milky Way, LMC, SMC, and other local galaxies. 

The reasonably bright gravitationally lensed quasar Q 0957 + 561 A,B at Zem — 1-41 
with separation 6" provides an excellent opportunity to study high-resolution spectra of two 
sightlines through a z > 1 DLA (e.g., Walsh, Carswell, & Weymann 1979; Wills & Wills 
1980; Young et al. 1980, 1981a,b; Turnshek & Bohhn 1993; Michalitsianos ct al. 1997; Zuo 
et al. 1997; Pettini et al. 1999). The z = 1.3911 absorber has a neutral hydrogen column 
density of A^(Hl) = 1.9 ±0.3 x 10^° cm'^ in the A spectrum and iV(Hl) = 8±2 x 10^^ cm'^ 
in the B spectrum (Zuo et al. 1997). This suggests that the DLA region is concentrated in 
front of the A image. At the redshift of the absorber, the physical separation of the A and B 
sightlines is 135h^^ pc (Smette et al. 1992). Zuo et al. (1997) found a differential reddening 
in spectrum A as compared to spectrum B and deduce a dust to gas ratio of ~ 0.6 in the 
DLA. Imaging of the quasar field has been reported by Dalhe, Maddox, & Lilje (1994, ground 
based) and Bernstein et al. (1997, WFPC2/HST). Despite these efforts, the luminous object 
associated with the z — 1.3911 absorber remains unidentified. 

The most extensive study of the DLA to date, using FOS/ HST spectra {R — 1300), 
was presented by Michalitsianos ct al. (1997). In general, they found "the differences in the 
damped Lyman series absorption in the lensed components are the only significant spectral 
characteristics that distinguishes the far-ultraviolet spectra of Q 0957 + 561 A,B." Assuming 
a spherical geometry and a neutral hydrogen fraction of A^(Hl)/A^(H) = 0.05 (Chartas et al. 
1995), the authors estimate the DLA mass to be 2.6 x 10^ Mq, which they interpreted to 
be consistent with a giant molecular cloud or a small condensation (spiral arm) within a 
galactic disk (viewed pole on). Because the metal lines are tightly correlated for the two 
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sightlines, Michalitsianos et al. (1997) suggested that the large metal-hne absorption arises 
in a galactic "halo" . 

In this paper, we present R = 45, 000 HIRES/Keck spectra of the Q 0957 + 561 A,B. 
We focus on the Mgll AA2796, 2803 doublet, the Mgl A2853 transition, and several Fell 
transitions arising in the DLA at z — 1.3911. In § 2, we present the data, and briefly 
describe the data reduction and analysis. We present the observed kinematic and spatial 
properties of the absorber column densities in § 3. Assuming photoionization equilibrium, 
we model the clouds and present their densities, sizes, and masses in § 4. We briefly discuss 
these properties in § 5 and provide concluding remarks in § 6. 

2. The Data 

The lensed quasar Q 0957 + 561 A.B was observed on the nights of 1994 November 
26-27 with the HIRES instrument (Vogt et al. 1994) on the Keck-1 telescope. HIRES 
was configured in first order using decker CI (7" in the spatial 0.861" in the dispersion 
directions, respectively). The resolution of R = 45,000 corresponds to a velocity resolution 
of FWHM = 6.6 km s^^. The resulting observed wavelength range is 4600 to 6900 A. As such, 
many different transitions were covered; however, we detected only the Mgll AA2796, 2803 
doublet, the Fell A2344, 2374, 2383, 2587, & 2600 multiplet, and the Mgl A2853 transition. 

The quasar images A and B were observed individually. Three 1500 second exposures 
were obtained for each image. The resulting spectra have signal-to-noise ratios of ~ 30 
per resolution element. The individual spectra were reduced and calibrated using the IRAF^ 
Apeoctract package for echelle data. The wavelength scale is vacuum and has been corrected to 
the heliocentric velocity. Continuum normalization was performed as described in Sembach 
& Savage (1992) and Churchill (1997). 

In Figure 1, we present the observed Mgll, Mgl, and Fen absorption profiles for the 
z — 1.391 absorber as a function of rest-frame velocity, where the A spectra (black) and B 
spectra (grey) are over plotted on the same velocity scale. The velocity zero point is set to 
z — 1.390861, which corresponds to the optical depth mean of the Mgll A2796 profile for the 
system in sighthne A (see Appendix A.l of Churchill & Vogt 2001). 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by AURA, 
Inc., under contract to the NSF. 
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3. Column Densities and Kinematics 

We modeled the data using Voigt profile (VP) decomposition. The fitting was performed 
with the code MINFIT (Churchill 1997; Churchill, Vogt, & Charlton 2002). The VP models 
provide the number of clouds and their velocities, column densities, and Doppler parameters. 
Because VP analysis is non-unique (especially in highly saturated profiles), and because it is 
sensitive to the signal-to-noise ratio and resolution of the spectra (Churchill 1997; Churchill, 
Vogt, & Charlton 2002), we use the VP parameters as a secondary means of studying the 
column densities and kinematics. Instead, we use the apparent optical depth method (Savage 
& Sembach 1991) as the primary means of studying the relationship between cloud velocities 
and column densities. 



3.1. Apparent Optical Depth Method 

Apparent column densities per unit velocity, Na{v) [atoms cm~^/(km s~-^)], were mea- 
sured for each transition using the formalism described by Savage & Sembach (1991). These 
column density spectra and their uncertainty spectra were then linearized to a common ve- 
locity binning of 2.23 km s~^ using fiux conservation. From these linearized data, an optimal 
column density was computed for each species in each velocity bin. The algorithm employed 
for computing optimal column densities has been described in Appendix A. 5 of Churchill & 
Vogt (2001). 

In brief, for adoption of the optimal Na{v) for an ion we employ one of three possibilities: 
(1) All transitions of an ion exhibit some saturation between Vb and so that Na{v) is a 
lower limit. Often, the transition with the smallest /A provides the best constraints on 
the lower limit; (2) All but one transition of an ion exhibits saturation, in which case, the 
adopted column density is taken from the unsaturated transition; and (3) All or more than 
one transition of an ion are unsaturated, providing multiple independent measurements. The 
optimal Na{v) is computed from the weighted mean in each velocity bin. 

The Mgll transitions in both systems exhibited unresolved saturation across the ma- 
jority of the profiles and therefore provide only upper limits on Na{v) at most velocities. 
Though the stronger Fell transitions (A2383 and A2600) exhibited saturation over the veloc- 
ity intervals (-30, +30), (+50, +70), and (+90, +110), the weakest transition (A2374), and 
to a lesser extent the A2587 transition, provide a robust Na{v) across these intervals. Mgl 
has no velocity regions with unresolved saturation. 

In Figure 2, we present N{v), the optimal apparent column densities, for Mgll (left). 
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Fell (center), and Mgl (right) as a function of rest-frame velocity^. The upper panels show 
N(v) for systems of sighthnes A and B (hereafter, systems A and B) as labeled and the lower 
panels show the column density difference, AN{v) = Na{v) —Nb{v). Limits are represented 
by arrows and 1 a uncertainties are shown as grey shading. 

A prominent feature of the data in Figure 2, is the N{v) peak for Fen and Mgl in system 
A at w ^ +20 km s~^. A second, even more prominent N{v) peak is at v ~ —10 km s~^ in 
system B. Interestingly, for Fell, the difference profile, AN{v), reveals positive (A>B) and 
negative (A<B) peaks with a quasi-pcriodicity of ~ 40 km s~^. In some cases these AA^(-;;) 
peaks arc due to different strengths of N{v) peaks aligned in velocity in both systems. In 
other cases the AN(v) peaks are due to the presence of an N{v) peak in the one system and 
the lack of an N{v) peak in the other at the same velocity. 



3.2. Velocity Correlations 

A striking feature of the data is the velocity alignment of an absence of absorption at 
V ~ +80 km s""*^, as can be clearly seen in the Mgll and strongest Fell profiles in Figure 1. 
This may indicate that the parcel of gas giving rise to absorption at f ~ +100 km s~^ is 
a separate physical entity from that giving rise to the lower velocity absorption. If so, this 
gas has httle to no velocity sheer across the sighthnes. At ~ —50 km s~^ and at ~ 
—100 km s"^ there are similar, yet less pronounced, profile inversions. In this case, there is a 
small difference of ~ 20 km s~^, which translates to a velocity sheer of 0.15/i75 km s~^ pc~^. 
Again, this could indicate that the higher (negative) velocity gas is physically distinct from 
the lower velocity gas. 

Ultimately, it is difficult to extract physical information directly from the flux values of 
such strong absorption lines, which exhibit saturation over much of the velocity interval. In 
order to further study the cloud-by-cloud kinematic connections between systems A and B, 
we ran a cross-correlation on the Fell and Mgl N{v) profiles. The cross-correlation function 
is defined by 

^ E [Na{v) - {Na)] ■ [Nb{v - Av) - {Nb)] 
a/E [Na{v) - {NA)f ■ E [Nsiv - Av) - {Nb)]' ' 

where Av is the lag velocity between the two systems. Equation 1 is defined so that a perfect 
correlation is ^ = +1 and no correlation is ^ = In Figure 3, we present ^ as a function of 



''Prom this point in the text, we drop the subscript "a" designating "apparent" in the AOD column 
densities. 



-8- 



Av. The left hand panels are Fell and the right hand panels are Mgl. The top (middle) 
panels show the self-correlation function for system A (B). By definition, these functions are 
symmetric about the lag velocity and have unity at zero lag velocity. 

The ^ for Fell in system A shows a remarkable pattern; there is a 40 km s~^ periodicity 
in the clouds, as evident in the two peaks in ^ at \Av\ ~ 40 and 80 km s~^. There is a lack 
of periodicity in the Fell for system B; the cross-correlation function reveals the velocity 
difference of the two strongest clouds separated by Av ~ 80 km s^"*^. The clouds in system 
A, in general, have large A^, so it is quite clear that this periodicity in system A is driving 
the shape of the AN{v) profile shown in Figure 2. 

The ^ for Mgl shows no periodicity, but only the velocity differences between the 
strongest components. Note the "noise" in C, for system B Mgl at 30 < |Af | < 60 km s~^, 
which has a magnitude of 5,^ ~ 0.1; this arises in the noisiest data and provides an estimate 
of the significance level of the stronger peaks for all of the cross-correlation functions. 

In the lower panels of Figure 3, we show the cross-correlation functions for system A 
against system B. The peaks in these functions provide the velocity difference between the 
strongest components in each system, which lies at Av ~ —27 km s~^ for both Fell and Mgl. 
A peak of ~ 0.7 is significantly below unity (approximately 3 a, based upon the above 
noise estimates) and quantifies the level at which the two Fell profiles do not resemble each 
other kinematically; ^ is dominated by the strongest components in complex profiles. 

Overall, this exercise reveals that there is no clear signal in the cross-correlation function 
for similar kinematics in the system A and system B profiles. This indicates that the clouds 
are not clearly traceable between the two sightlines. 

3.3. Integrated Column Densities 

In Table 1, we present the integrated apparent column densities, A^, for the Mgll, Fell 
and Mgl ions. The A^ were computed for fixed velocity intervals (from Vi to V2) using the 
data presented in Figure 2. The velocity intervals were defined by local N{v) minima in the 
Fell spectra for system A and B individually. Six velocity intervals were found for system 
A and seven were found for system B. These intervals roughly represent individual "clouds" 
giving rise to the complex absorption profiles. 

As stated above, the Mgll transitions provide only upper hmits for most velocity inter- 
vals. The exception is velocity interval 1 (or "cloud number 1") in system A. Clouds 2, 3, 5, 
and 7 in systems B are marginally saturated; there is at least one saturated pixel in each of 
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these clouds. For this reason, we quote these particular clouds in system B as upper limits. 
However, unless there are high column density clouds with b < 2 km s~^ at the location of 
these saturated pixels, the quoted values for system B could be marginally acceptable as 
measurements. We choose to not invoke them as measurements for our analysis. 

For both systems, we computed the column density ratios A^(Mgl)/A'"(Fell) for each 
cloud and listed them in Tabic 1. The A^(Mgl)/A^(Fell) data are plotted in Figure 4. The 
horizontal bars of the data points give the velocity interval of the N{v) integrations and the 
vertical bars give the uncertainties in the column density ratio. Plotted in the upper panels 
of Figure 4, and aligned in velocity space for ease of inspection, are the Fell and Mgl N{v) 
profiles. 

As seen in Figure 4, the ratios for all clouds in both systems A and B are consistent 
with logiV(Mgl)/iV(FeIl) ~ —2. This result suggests a high level of uniformity in both the 
gas-phase [Mg/Fe] abundance patterns^ and ionization conditions across velocity space. 

The VP subcomponents describing these dominant "clouds" also have ratios relatively 
consistent with these findings, with —2.02 ± 0.04 in system A and — 2.1 ± 0.3 in system B. 

Integration of N{v) across the full profiles yields log A^(Mgl)/Ar(Fell) = -2.03lo:o8 for 

system A and — 1.961°' {5 for system B. These quantities arc shown as grey shaded regions on 
Figure 4. The individual clouds have ratios that are consistent with that of the full system. 
This level of uniformity with velocity is consistent with that reported by Prochaska (2002) 
for a sample of 13 higher rcdshift DLAs. Worth noting is that these data also provide direct 
evidence of spatial uniformity on the scale of 135h^^ pc in a DLA. 

4. Nature of the Absorbers 

4.1. Photoionization Modeling 

Inferring the physical conditions in the clouds is model dependent. Here, we assume 
photoionization equilibrium. Using the measured column densities as constraints on the 
models, we derive the cloud metallicities, densities, sizes, and masses. 

We used Cloudy (Ferland 1998) to construct grids of model clouds. We assume that the 
A^(Hl) in system A clouds is logiV ~ 20.3 cm~^ and in system B clouds is log A?" = 19.9 cm~^ 
(Zuo et al. 1997). Since the measured A^(Hl) is actually the sum of individual clouds in a 



'We use the notation [X/F] = log{X/Y) — log(X/y)0 throughout this paper. 
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system, we note that some inferred properties (sizes and masses) for the individual clouds 
will be overestimated. 

Each cloud is modeled as a constant density plane-parallel "slab" with ionizing radiation 
incident on one face. For the ionizing flux, we used the z — lA ultraviolet background (UVB) 
spectrum of Haardt & Madau (1996) and Madau, Haardt, &: Rees (1999). The separation 
velocity of the DLA and the quasar is ~ 2500 km s~^. Thus, it could be argued that the 
absorber is not an intervening system, but is associated with the quasar itself. If so, the 
quasar flux, and not the UVB flux, would dominate the photoionization of the gas. 

One indicator of associated absorption is the presence of partial covering (Barlow & 
Sargent 1997; Hamman 1997; Ganguly et al. 1999). None of the low ionization transitions 
we observed show signs of partial covering. That is to say, that the fully saturated regions 
of the Mgll AA2796, 2803 doublet members are black (statistically consistent with zero flux) 
in their cores. Also, at the same velocity intervals where Fell A2383 and A2600 are also fully 
saturated, their absorption is also black. As such, we assume that the DLA is intervening. 
If the absorber is intervening, then the cosmological separation (Az = 0.02) from the quasar 
is ~ IS/ifg^ Mpc. It has been argued that, at this distance, the quasar flux does not affect 
the ionization balance in the absorber (Michahtsianos et al. 1997). 

At hv — 13 eV, the flux normahzation is uFi, = 1.26 x 10~^ erg s~^ cm~^. We output 
the column densities for selected ionic species (esp. Mgll, Mgl, and Fell) as a function 
of ionization parameter, U — n-y/rin. This quantity is the ratio of the number density of 
photons capable of ionizing hydrogen to the number density of hydrogen atoms. We explored 
the range — 8 < logC/ < —1 in 0.1 dex intervals. 

For these large A'"(Hl) values, the model clouds are optically thick and have an extended 
neutral layer. As such, the relative dependence of the column densities with ionization 
parameter are indistinguishable between the logN — 20.3 cm~^ and logN — 19.9 cm~^ 
cloud models. 



4.2. Metallicity 

In Table 2, we hst the 3 a rest-frame equivalent widths of selected transitions. We infer 
upper limits on the cloud metallicities using the Znll A2026 transition, which only weakly 
depletes onto dust and is known to trace Fe-group elements (e.g.. Savage & Sembach 1996; 
Lauroesch et al. 1996). No Znll was detected in our spectra; the 3 a equivalent width limit 
was 13 mA in system A and 15 mA in system B. Without detections of CrII and Znll we 
cannot directly estimate the effects of dust in these systems (see Zuo et al. 1997). 
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For these equivalent widths, Znll A2026 is effectively on the linear portion of the 

curve of growth. We obtain an upper limit of logA^(Znll) < 11.9 cm~^ for system A and 
logA^(Znll) < 12.0 cm~^. In the cloud models iV(Znll) is affectively independent of ioniza- 
tion parameter for log?7 < —3.5 and decreases by only 0.5 dex from logU = —3.5 to — 1. 
For logf/ < —3.5, we obtain an upper limit of [Zn/H]< —1 for system A and [Zn/H]< —0.5 
for system B. 

These values are corroborated by upper limits based upon the weaker Znll A2063 tran- 
sition. Our HIRES/Keck spectra provide only slight improvements over the limits (also 3 a) 
of [Zn/H]< —0.75 for system A and [Zn/H]< —0.31 for system B reported by Pettini et al. 
(1999). 

4.3. Ionization Parameter: A Matter of Dust 

In the model clouds, A^(Mgll) and A^(FcIl) are virtually independent of U, whereas 
iV(Mgl) decreases about 0.8 dex for every 1 dex increase in U. The A^(Mgl)/A^(MgIl) ratio 
could provide an estimate of the cloud ionization parameters; however, the Mgll profiles 
are severely saturated. Since the N{v) profiles for Fell are well defined at all velocities, the 
A^(Mgl)/A^(Fell) ratio, in principle, holds potential for inferring the ionization parameter. 
In Figure 5, we show the Fell and Mgl column densities as a function of logC/ for one of 
our Cloudy models. This model cloud has logA^(Hl) = 20.3 cm~^ and a solar abundance 
pattern. For presentation purposed, we have scaled the metallicity to [Z/Zq] = —2 to 
match the integrated N for system A presented in Table 1. The logarithm of the ratio 
iV(Mgl)/A^(FeIl) decreases from ^ —0.5 at logf/ = — 8 to ~ —4.5 at \ogU = —3, after 
which it begins to decrease rapidly. For a solar abundance pattern and no dust depletion, 
the cloud models for both systems have log A^(Mgl)/A^(FeIl) ~ —2 for logf/ ~ —6. 

The fact of the matter, however, is that A^(Mgl)/A^(FeIl) is directly proportional to the 
gas-phase [Mg/Fe] abundance ratio, which can be strongly affected by the dust depletion 
factors for Mg and for Fe. The study of the effects of dust in DLAs is one of central 
importance (e.g., Ellison et al. 2001, and references therein). 

Zuo et al. (1997) report a dust to gas ratio of ~ 0.6 for the z = 1.3911 DLA studied here. 
However, we note that there are several assumptions invoked for their calculations, including 
extrapolation of the Galactic extinction law and no time delay or variability between the 
lensed images. The latter is certainly not strictly true {At ~ 420 days, Kundic et al. 1997; 
Haarsma et al. 1999), and the former may not apply to DLAs (e.g., Pettini et al. 1997; 
Prochaska 2002) even if the extrapolation accurately represents the Galactic law. Even 
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without a definitive estimate of the dust content in this DLA, we can explore how dust 
depletion might effect the inferred cloud properties. 

In the Galaxy, dust depletion factors for Mg range from S — —0.3 in warm "halo" gas to 

5 — —1.6 in cool "disk" gas and for Fe range from 5 — —0.6 to —2.3 (Lauroesch et al. 1996; 
Savage & Sembach 1996; Welty et al. 1999, 2001). In order to allow for the observed range 
in [Mg/Fc] due to dust depletion, we also consider cloud models with warm halo and cool 
disk dust depletion factors. These ranges also encompass the observed range of a-group to 
Fe-group variations observed in the photosphere of Galactic stars (McWilliam et al. 1995; 
Lauroesch ct al. 1996; Johnson 2002), globular cluster stars (Smith ct al. 2000; Stephens 

6 Bosegaard 2002), LMG/SMC stars (Venn 1998), and dwarf galaxies (Shetrone, Cote, & 
Sargent 2001). 

For the warm halo pattern we find the lower limit of log U ~ —6.2 and for the cold disk 
pattern we find the upper limit of logU ~ —5.1; we have. 



for both the system A and system B clouds, where the no-dust, solar abundance pattern 
value is log U — —6.0. 

We adopt this range in logC/ as an estimate of the inferred ionization condition in 
the clouds. We note, however, there is growing observational evidence that dust depletion 
in DLAs may not be significant and may be fairly uniform both from system to system 
(Pettini et al. 1999; Ellison et al. 2001) and from velocity component to velocity component 
within a system (Prochaska 2002). However, we note that there are examples of possible 
intrinsic abundance variations of [Mn/Fe] in DLAs (Ledoux, Bergeron, & Pctitjean 2002) 
and of strong dust depletion in DLAs selected because they are strong molecular hydrogen 
absorbers (e.g., Petitjean, Srianand, & Ledoux 2002). 

The statistical data supporting low dust content in DLAs would suggest that the above 
range of ionization parameters serves as a somewhat conservative approach for bracketing 
the inferred ionization parameter. The largest uncertainty lurks in the original assumption 
that the clouds are in photoionization equilibrium. 



For the UVB spectrum normalized ai z = 1.4, the relationship between the hydrogen 
number density of the clouds and the ionization parameter is 



-6.2 <\ogU < -5.1 



(2) 



4.4. Densities, Sizes, and Masses 



log riH = —4.9 — log U cm 



-3 



(3) 
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for constant density cloud models. For the inferred range of ionization parameters, we find 

3 < nn.A < 20 cm'^ and 2 < nn,B < 15 cm^^, (4) 

for the system A clouds and system B clouds, respectively. These densities fall within the 
lower density range for cold clouds in the Galactic interstellar medium (Spitzer 1985; Savage 
& Sembach 1996). The cloud line-of-sight physical extent, D — N{}i)/nu, is then, 

log D = log iV(H) + log ^7- 13.7 pc, (5) 

where A'"(H) is the total hydrogen column density. For the inferred range of ionization con- 
ditions, the ionization fraction of hydrogen is negligible so that the approximation A'"(H) = 
N{m) holds^. We obtain, 

5<Da<25 pc, and 1 < < 15 pc, (6) 

for the system A and system B clouds, respectively. To the extent that the photoionization 
modeling has provided a reasonable approximation of the cloud physical conditions, we can 
infer that the cloud sizes are roughly a factor of ten smaller than the sightline separation of 
the lensed quasar. 

We can estimate the cloud mass under the assumption of spherical symmetry, which 
gives M (X miinii{D/2Y, where mn is the mass of hydrogen. In terms of our parameteriza- 
tions, the cloud mass is, 

logM = 31ogA^H + 21ogC/-47.7 Mq, (7) 
which alternatively can be written as 

Mc^2Nl,U'_, Mo, (8) 

where A^2o is the total hydrogen column density of the cloud in units of 10^° cm~^, and C/_6 
is the ionization parameter in units of 10~^. We find approximate masses of 

10 < Ma/Mq < 1000 and 0.5 < Mb/Mq < 60, (9) 

for the system A and system B clouds, respectively, for the range of ionization parameters 
given in Equation 4. 



^The adopted Cloudy models yield hydrogen ionization fractions slightly smaller than the value 0.05 
reported by Chartas et al. (1995). 
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If we instead assume that the clouds are cyhndrical "slabs" with "height" D — N{E.)/nH 
and "radius" R, we have M oc munuDR'^, which can be simphfied to 

Mc^25N2oRloo Mo, (10) 

where i?ioo is the cylinder radius in units of 100 pc. Note that the mass for this geometry 
is independent of the cloud density, nn, and therefore ionization parameter, U. Given that 
we cannot track the individual clouds from sightline A to sightline B, a conservative upper 
limit of Rioo — 2hj^ can be applied. With this upper limit, we obtain 

Ma < mh^^ Mq and Mb < 160h^^ Mq, (11) 

for upper limits on the system A and system B cloud masses. 

4.5. Caveats 

We examined the effects on the above inferred cloud properties of different shapes of 
the ionizing spectrum. We examined a grid of Cloudy models in which O stars and B stars 
contributed to and/or dominated over the UVB. We followed the formalism of Churchill, 

Vogt, & Charhon (2002) and Churchill & Le Brun (1998) (esp., see Figures 12 and 13 of 
Churchill, Vogt, & Charlton (2002)). We find that the contribution of O and B stars makes 
virtually no difference in the inferred cloud properties. This is due to the high level of self- 
shielding in the clouds. Thus, we find that the inferred cloud properties are robust under 
the assumption of photoionization equilibrium. 

It is important to point out that the above inferred sizes and masses are overestimates 
within the formalism we have utilized. Each individual cloud was assumed to have identical 
A^(Hl), and the value applied was the total A^(Hl) for each respective system. That is, all 
system A clouds were assumed to have log A^(Hl) ~ 20.3 cm~^ and all system B clouds were 
assumed to have log A'"(Hl) = 19.9 cm~^. Therefore, the inferred sizes and masses for these 
clouds are smaller than we have quoted here. 

5. Discussion 

As first inferred from the FOS/ HST spectra and confirmed here, the absorber appears 
to have a projected transverse physical extent no greater than ~ 135h^^ pc and seems to 
be enshrouded by an optically thick gaseous complex (Michalitsianos et al. 1997). We have 
resolved the kinematics of the low ionization gas in these systems and find the total velocity 
spread to be ~ 300 km s'^, based upon the Mgll AA2796, 2803 doublet. 
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We have defined six "clouds" for system A and seven clouds for system B. The data and 
photoionization models suggest a picture in which the sightline physical extent (sizes) of the 
clouds in this DLA are of the order 10 pc. The inferred densities are ~ 10 cm~^, and the 
masses are no greater than ~ AOOh^^ Mq. The temperatures of the cloud models, if they are 
taken at face value, have a gradient with values that range from T ^ 2000 K at the cloud 
face to T ~ 500 K at the cloud core. To the extent that the assumption of photoionization 
equilibrium is appropriate and has been modeled accurately, we have found that the DLA 
clouds are similar to low mass (10-100 Mq), low density (1-10 cm~^), cold (~ 1000 K) Hi 
clouds in the Milky Way (Spitzer 1985; Savage & Sembach 1996). 

5.1. Where lies the DLA? 

Directing our attention to the physical nature of DLAs, we focus on two questions. (1) 
Is the DLA region itself only one of these small clouds, and therefore does it exhibit very 
narrow kinematics? (2) Is the DLA region only a single phase of gas, or can it arise in a cold 
phase and a warm phase (e.g.. Lane, Briggs, & Smette 2000)? The answer to these questions 
will provide clues for interpreting the observed abundance ratio and ionization uniformity 
between the clouds. 

Based upon the cross-correlation (see Equation 1) of the system A and system B pro- 
files, we find that the individual Mgll and Fell clouds are not identifiable between the two 
sightlines separated by 135hj^ pc. This implies that we do not track the same clouds from 
sightline A to B, which in turn implies that they are smaller than the line of sight separa- 
tion. Ranch et al. (2002b) found a similar lack of a clear identification between clouds in an 
optically thick Mgll absorber at z = 0.5656 in the triple sightlines of Q 2237 -|- 0305, which 
have separations of 135-200/i75^ pc. 

However, the cross-correlation function has a single strong peak, indicating that there is 
a single strongest component in each system; system A has a strongest absorbing component 
at f ~ +20 km s~^, and system B has a dominant component at f ~ —10 km s~^ (see 
Figure 2). It is reasonable to argue that the DLA region is physically associated with 
the sites of strongest absorption. In numerical simulations of hierarchical clustering, the 
strongest absorption lines arise in the single most highly peaked baryon overdensities (e.g. 
Haehnelt, Steinmetz, & Ranch 1996; Ranch, Haehnelt, & Steinmetz 1997). 

Because there are only low resolution spectra available for the neutral hydrogen lines 
(e.g. Turnshek & Bohlin 1993; Michalitsianos et al. 1997; Zuo et al. 1997; Pettini et al. 
1999), the measured A'"(Hl) contains no direct kinematic information. However, using the 
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kinematics of the Fell gas as a template for the neutral hydrogen, we use the Lyct, Ly/3, Ly7, 
LyS, and Lye profiles from the FOS spectrum in the HST archive to test the hypothesis 
that the Hi arises in a single cloud. While the Lya line constrains the damping wings, the 
higher series lines provide constraints on the Doppler parameters. 

We ran three simulations on the system A profiles: (1) a single cloud at the redshift of 
the strongest Fell-Mgl component; (2) six clouds with equal A^(Hl) with velocities aligned 
with the Fell-Mgl components; and (3) six clouds with iV(Hl) proportional to A^(Fell) and 
velocities aligned with the Fell-Mgl components. For simulations "2" and "3", we assumed 
all clouds had the same Doppler parameter. The single cloud model fit the data well; because 
of the higher order Lyman series, we were able to put a limit of 50 < 6(Hl) < 60 km s~^ 
on this single cloud. Although both multiple cloud simulations yielded good fits to the 
higher order Lyman series, they totally failed to fit the Lya damping wings associated with 
logA^(Hl) =20.3 cm-2. 

Using the simulations as a guide, it is reasonable to assume that the measured A'"(Hl) 
for each sightline is the value associated with the strongest component in systems A and 
B. One scenario is that these strongest clouds are physically hnked. If so, then the velocity 
shear in the DLA is Av/Ar ~ 30/135 = 0. 22/175 km s^^ pc~^ and (for a planar geometry) 
the spatial column density gradient is ~ 6 x 10^^ cm~^ pc~^. It is more likely, however, that 
the clouds arc individual parcels of gas. If so, then the A and B sightlines probe physically 
distinct clouds with different A'"(Hl) reflecting the structure of the overall absorbing complex. 



5.2. A Two-Phase DLA? 

Using 21-cm observations of the z — 0.0912 system toward the radio quasar B 0738+313, 

Lane, Briggs, & Smette (2000) reported a two phase DLA with Tcoid — 300 K and Twarm — 
5000 K, where the warm phase contributes roughly two thirds of the A^(IIl). In the case of 
Q 0957 + 561, VLBI mapping of the quasar (e.g., Garrett ct al. 1994; Campbell ct al. 1995) 
reveals that the extended radio emission is probably not covered by the DLA. Kanckar & 
Chengalur (2002) reported no detection of 21-cm absorption with an RMS of 2.4 mJy at a 
resolution of 3.9 km s~^. These facts are consistent with the DLA having a small covering 
factor. It may not be possible to use 21-cm absorption to constrain the phase structure of 
the DLA toward Q 0957 + 561. 

However, based upon the VP fitting to the dominant components of the Mgl profiles 
of the z — 1.3911 DLA, and as can be seen in Figures 2 and 4, the Mgl profiles have 
broad wings. The VP fits of these features yielded a single broad component with b — 
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15.3 ± 1.7 km centered at = 17.2 km in system A with logA^(Fell) = 14.10 ± 
0.02 cm~^ and logA^(Mgl) = 12.08 ± 0.05 cm~^. In system B, three subcomponents were 
found. A narrow component with b = 5.6 ± 1.8 km s^^ centered at w = —3.3 km s~^ with 
log A^(Fell) = 13.92 ± 0.04 cm^^ and log A^(Mgl) = 11.7 ± 0.1 cm~^, a broader component 
with 6 = 10.2±3.9km s'^ centered atv = -14.9 km s^^ with log iV(Fell) = 13.73±0.02 cm'^ 
and logiV(Mgl) = 11.7±0.2 cm~^, and a third broader component with b = 7.6±0.4 km s~^ 
centered at ^; = +14.0 km s~^ with logiV(Fell) = 13.40 ± 0.02 cm^^ and logA^(Mgl) = 
11.3 ±0.4 cm-2. 

The broadness of the VP component in system A and the two adjacent VP subcompo- 
nents in the wings of the system B component may be revealing a warm phase analogous to 
that reported by Lane, Briggs, & Smette (2000). If such ionization structure were present, it 
would impact the inferences we have on the cloud sizes and masses from our photoionization 
models. We ran a fourth simulation on the Lyman series lines in which we model a "cold" 
phase cloud and a "warm" phase cloud, each contributing half the total A^(Hl) at the redshift 
of the strongest Fell-Mgl component. 

The maximum Doppler parameter the warm phase can have is 6 = 60 km s"^ as con- 
strained by the higher order Lyman series lines. We adopted 6 = 55 km s~^. By itself, this 
"warm" component, with logiV(Hl) = 20.0 cm~^, cannot generate the Lya damping wings 
and also does not fit the Ljf3 wings. Testing a range of Doppler parameters in the range 
5 < 6 < 15 km s^^ for the "cold" phase, we find that this component alone cannot account 
for a substantial portion of the absorption in any of the Lyman scries lines. However, to- 
gether the two components provide a fit to the data that is statistically identical to the single 
component model, with 6 < 10 km s"-*^ in the "cold" component yielding the best fit. 

Therefore, we can confidently state that the Hi arises at a single velocity, or velocity com- 
ponent; however, we cannot distinguish between a single phase model and a "cold+warm" 
two-phase model. We note however, that the inferred temperatures of our two-phase 
model are higher than the temperatures reported by Lane, Briggs, & Smette (2000) for 
the z = 0.0912 system toward B 0738 + 313. 

It is difficult to incorporate a two-phase model into the direct observation that the 
ion ratios are exceptionally uniform from cloud-to-cloud across the full velocity extent of 
the profiles. We also find that this uniformity is spatial on the scale of ISbh^^ pc. As 
stated above, this would suggest that the gas-phase abundances of [Mg/Fe] are uniform 
both kinematically and spatially vis-a-vis the ionization conditions. 

Observational evidence is that Lyman limit systems are multiphase, and they do not 
show kinematic uniformity. Examination of A'"(Mgl)/A'^(Fell) ratios in the Lyman limit 
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systems studied by Churchill & Vogt (2001) reveals significant cloud to cloud variation with 
velocity in the integrated apparent column densities (see their Table 6). For example, the 
z — 1.0479 Mgll system in PG 0117 + 213 has ratios —1.03 and —1.55 in its two kinematic 
subsystems. The z = 0.7729 system in PG 1248 + 401 has a remarkable variation from —1.92 
to —0.56 in its two subsystems. In addition, there are examples of strong abundance and/or 
ionization variation with kinematics (e.g., Ganguly, Churchill, & Charlton 1998) and with 
spatial separation (e.g., Rauch et al. 2002b) in sub-Lyman limit systems. 

5.3. Systematic Kinematics? 

The kinematics of the Mgl profiles are reminiscent of the asymmetric profiles studied 
by Prochaska & Wolfe (1997, 1998), who promoted the idea that DLA gas arises in thick 
rotating disks of galaxies. However, in the case of the z = 1.3911 DLA toward Q 0957 + 561 
the inferred size of the individual cloud giving rise to this absorption is ~ 10 pc and this 
precludes that this particular profile asymmetry is generated by disk kinematics. 

It is noteworthy that the velocity sheer of the strong absorbing components is similar 
to the velocity sheer of the profile inversions (regions where there is a paucity of gas; sec 
first paragraph in § 3.2). This fact is indeed suggestive of some coherence in the overall 
absorbing structure and dynamics, even if the individual clouds cannot be directly traced 
across sightlines. The profile inversions may indicate physically separated gas parcels, as 
suggested by the small cloud sizes we have derived. If so, then the common velocity sheer 
is remarkable; the multiple absorbing complexes would be members of a generally extended, 
kinematically systematic object (a co-rotating halo?; Weisheit 1978; Steidel et al. 2002). 

It is clear from the higher ionization data (Michalitsianos et al. 1997), that there is at 
least a low density, high ionization phase associated with this system in addition to that 
giving rise to the Mgll, Fell, and Mgl. It is likely that this gas, especially the Civ, is more 
diffuse and extends more smoothly across the sightlines (Rauch, Sargent, & Barlow 2001), 
and is therefore not as directly coupled to the low ionization systematic kinematics. 

We examined if the velocity difference of ~ 30 km s~^ between the strongest Mgl clouds 
in the A and B profiles could arise due to disk rotation. Using the formalism of Charlton, 
Churchill, & hinder (1995), we modeled the velocity sheer for a ~ 200 pc separation scale due 
to disk kinematics for random sightline orientations. We assumed a disk circular velocities 
in the range of 150 < Vc < 250 km s~^. We find that the sightlines cannot have velocity 
differences as large as 30 km s~^ but for highly contrived sightline-galaxy orientations (the 
sightlines are too close together) . 
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However, our models do not incorporate gas dispersion perpendicular to the plane of 
the disk (e.g., Charlton & Churchill 1998). Since such motion is most definitely present in 
galactic disks, these models by no means rule out the possibility that the absorbing material 
arises in a disk; they do strongly suggest that the velocity sheer between the two sightlines 
is not a result of simple disk kinematics. 

6. Conclusion 

Wc have observed the images of the lensed quasar Q 0957+561 A,B with the HIRES/Keck- 
1 instrument (resolution FWHM ~ 6.6 km s~^). Wc have presented an analysis of the Mgll, 
Mgl, and Fell absorption profiles from the z = 1.3911 DLA system. We adopted the hydro- 
gen column densities for systems A and B from Zuo ct al. (1997), which are log A^(Hl) = 20.3 
and 19.9 cm^^, respectively. The line of sight separation is ~ 135h^^ pc at the redshift of 
the absorber (Smette et al. 1992). 

We converted the absorption profiles to their apparent optical depth column density 
{N{v), Savage & Sembach 1991). Based upon the location of local minima in the N{v) 
profiles for Fell, we defined six "clouds" in system A and seven clouds in system B and 
integrated the N{v) to obtain the cloud column densities. There is a "dominant" cloud in 
each line of sight. It may be that these clouds contain the bulk of the neutral hydrogen gas. If 
the cloud geometry is planar and extents across sightline B, then there is a neutral hydrogen 
column density gradient of 9 x 10^^/175 cm~^ pc~^ and a velocity sheer of ~ 0.2/175 km pc~^. 

The clouds were assumed to be in photoionization equilibrium. Using Cloudy (Ferland 
1998), we modeled the clouds as constant density, plane-parallel "slabs" illuminated on 
one face by the ultraviolet background ionizing spectrum. We used the A'"(Mgl)/A'^(Fell) 
ratio in each cloud to constrain the ionization conditions. Since both Mg and Fe suffer 
dust depletion and originate predominantly in separate nucleosynthctic environments, we 
bracketed the [Mg/Fe] abundance pattern for the range of dust depletions seen in the Galaxy 
and LMC/SMC and for the observed abundance patterns in the local universe. 

The observed A^(Mgl)/iV(Fell) ratio is remarkably uniform. Not only are the ratios 
consistent with log[A^(Mgl)/A^(FeIl)] = —2 across the full velocity range in both systems, 
but they are also consistent with this value spatially (in both sightlines). This resulted in 
very uniform cloud physical properties as inferred from the photoionization modeling. The 
ionization parameter of the clouds is in the range —6.2 < logU < —5.1. This yielded clouds 
with densities of 2 < nn < 20 cm~^ and fine of sight physical extents of 1 < D < 25 pc. The 
inferred masses are geometry dependent. For spherical geometries the masses of the clouds 
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in system A are 10 < M/Mq < 1000 and in system B are 1 < M/Mq < 60. For cylindrical 
geometries constrained by the line-of-sight separation of less than 200 pc, the cloud masses 
have upper limits of AQQh^^ and IGO/ifg^ for systems A and B, respectively. These 
cloud properties are consistent with those for lower density, cold clouds in the Galactic 
interstellar medium (Spitzer 1985; Savage & Sembach 1996). 

We focused our discussion on the physical nature of the DLA "region", the object 
that actually gives rise to the damped Ly« absorption of logA^(Hl) = 20.3 cm^^. Based 
upon simulations, we favor a picture in which the DLA is a single cloud in the multi-cloud 
profiles. We cannot discern, however, if the DLA comprises a "cold" single ionization phase, 
as suggested by our photoionization models, or a "cold+warm" two-phase gas complex. 

If the DLA cloud is spherical in nature, then its size is on the order of ~ 10 pc, and 
it is limited to one of the sightlines (A). This implies a covering factor of less than 0.1. 
The other multiple gas clouds in the proximity of this small DLA cloud would have to have 
experienced the same sources of nucleosynthetic enrichment, be optically thick in A^'(Hl), 
and have similar dust contents. This implies that the material distributed in proximity to 
the DLA is well mixed and ionized uniformly. This is in stark contrast to the significant 
variations seen in Lyman limit systems (e.g., Churchill & Vogt 2001), which are thought to 
arise in the outer disks and halos of galaxies. As such, we suggest that the low ionization 
clouds accompanying DLAs are not arising in galactic halos. 

Rather, we infer that DLAs arise in small gas-rich regions within galaxies. The data and 
models suggest that these regions are complexes comprised of small, optically thick clouds 
similar to the lowest mass, cold Hi clouds in the Galaxy. Furthermore, the data suggest that 
they are well mixed chemically and have similar photoionization conditions. 

Support for this work was provided by NASA (NAG 5-6399), and by the NSF (AST 
96-17185). We thank Mike Keane for his assistance with the observations. We also thank 
Patrick Hall, Don Schneider, and an anonymous referee for insightful comments that helped 
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Table 1. Apparent Column Densities^ 



System A 



Cld# 


1 


2 


3 


4 


5 


6 


{VuV2) 


(-140,-88) 


(-88, -42) 


(-42, -3) 


(-3,35) 


(35,80) 


(80, 124) 


iV(Mgll) 


. o 24+0.02 


> 13.5 


> 13.6 


> 13.5 


> 13.6 


> 13.5 


iV(Fell) 


1 q 9Q+0.02 
iO.ZO_o.o2 


1 q p: 7+0.02 
lO.O ( _o.02 


13.701°;^^ 


io qr+0.03 


io f^q+0.02 


1 s 64+0-03 


iV(Mgl) 


< 10.6 


11 fin+0.09 


11.541H? 


1 2 04+°-°^ 

iZ.U^_0.04 


11 arr+0.08 
ii.UO_Q 09 


11.56+H? 


Mgl/Fell 


< -2.6 


_i 07+0.09 
^■^ ' -0.10 


9 1^+0.09 
^•-'-O-O.ll 


_i q 1+0.05 

-'-•^-'--0.04 


^••^^-0.09 


-2.08in? 



System B 



Cld# 


1 


2 


3 


4 


5 


6 


7 


(^^1,1^2) 


(-140, -90) 


(-90,56) 


(-48,-16) 


(-16, +5) 


(+5, +40) 


(+40, +85) 


(+85, +125) 


Ar(Mgll) 


> 13.5 


> 13.4 


> 13.4 


> 13.3 


> 13.4 


> 13.6 


> 13.3 


N{Feu) 


13.471°:°^ 


1 97+0-02 
±0.z / _o.o2 


13.40iHi 


1 88+0-05 


13.48lH^ 


13.481H^ 


13.07tH^ 


N{Mgi) 




11 28+°-^^ 

ii.ZO_o 20 


ii.47tS:}^ 


11 87+°-°^ 

-'--'-•O' -0.05 


11-651°:°^ 


11.641°:°^ 


ii-3ilg:i^ 


Mgl/Fell 


9 qo+0.22 


_1 qq+O.1.5 
-'-■■-'•-'-0.20 


-i.93;^-|^ 


_9 r)i+0.07 


-1.83;^;[5« 


1 Q 1+0.09 


-i.76;[5:- 



All column densities are in units of atoms cm ^. 
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Table 2. Equivalent Width Limits of Selected Transitions'^ 



Transition LOS A LOS B 

Wr,lim A Wr,lim A 



Mgl 2026 


0. 


,012 


0. 


,016 


Znll 2026 


0. 


,013 


0. 


,015 


Cm 2056 


0. 


,014 


0. 


,016 


Cm 2062 


0. 


,011 


0. 


,014 


Znll 2063 


0. 


,011 


0, 


,013 


Cm 2066 


0. 


,010 


0. 


,012 


Fell 2250 


0. 


Oil 


0. 


,013 


Fell 2260 


0. 


,012 


0. 


,013 


Nil 2312 


0. 


,010 


0. 


,013 


Nil 2321 


0. 


,012 


0. 


,014 


Fel 2484 


0. 


,010 


0. 


,013 


Sil 2515 


0, 


,012 


0, 


,015 


Fel 2524 


0. 


,010 


0. 


,014 


Mnll 2577 


0. 


,010 


0. 


,013 


Mnll 2594 


0. 


,027 


0. 


,017 


Mnll 2606 


0. 


Oil 


0. 


,014 



^AU limits are rest-frame 3 a. 
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Fig. 1. — Normalized HIRES/Kcck spectra of the Mgll, Mgl, and Fell absorption profiles 
presented in the system rest-frame velocity. A spectra are black and B spectra are grey. 
The velocity zero point is set to 2; = 1.390861, which is the optical depth mean of the Mgll 
A2796 transition in the A spectrum. 
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Fig. 2. — The optimal apparent column densities profiles, N{v), for Mgll (left), Fell (center), 
and Mgl (right) presented in the system rest-frame velocity. The grey shading is the 1 a 
uncertainties in the N{v) profiles. The lower panels show the differences A— B. 
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Fig. 3. — The cross-correlation function, ^(Af), for the Fell and Mgl N{v) profiles. Left 
hand panels are Fell with the self-correlation of A (top), self-correlation of B (middle), and 
cross-correlation of A with B (lower). The right hand panels are Mgl. 
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Fig. 4. — The optimal apparent column densities profiles, N{v), for Fell and Mgl presented 
in the system rest-frame velocity. A spectra are on the left and B spectra on the right. The 
bottom panels show the ratios of the integrated column densities for each of the velocity 
bins corresponding to "separate clouds" in the Fell profiles. The dotted horizontal line is 
the average of the individual clouds and the grey shaded area is the standard deviation. This 
average and standard deviation are written in the corner of the respective panels. 
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Fig. 5. — A photoionization equilibrium Cloudy model of a plane-parallel cloud with flux 
illuminated on one face. The cloud has logiV(Hl) = 20.3 cm~^, a metallicity of [Z/Zq] = —2 
and a solar abundance pattern (i.e., no dust depletion). Shown are A^(Fell) and A^(Mgl) as 
a function of ionization parameter (see text) . 



